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For  round  shaped  heat  source  or  heat  sink,  the  annular  thermoelectric  couple  (ATEC)  with  annular  shaped 
legs  instead  of  flat-plate  thermoelectric  couple  (FTEC)  is  proposed  to  solve  the  performance  deterioration 
caused  by  geometries  mismatching.  Based  on  one-dimensional  steady  model,  the  fundamental  formulas 
of  ATEC  are  derived.  The  impact  of  geometric  feature  of  ATEC  manifested  by  the  annular  shaped  param¬ 
eter  sr,  on  the  dimensionless  output  power  and  efficiency  of  ATEC  under  different  temperature  ratios, 
external  loads  and  dimensionless  figure  of  merits  is  examined  theoretically.  Also,  to  enhance  the  practi¬ 
cability  of  ATEC,  some  suggestions  or  notes  are  provided.  Results  show  that  the  forms  of  fundamental  for¬ 
mulas  about  ATEC  are  the  same  as  those  about  FTEC.  Compared  to  temperature  ratio  and  dimensionless 
figure  of  merit,  the  influence  of  sr  is  affected  by  external  load  seriously.  The  extremums  of  dimensionless 
output  power  and  efficiency  are  presented  on  condition  that  sr  =  1.  As  the  closer  sr  is  to  1,  the  larger  the 
dimensionless  output  power  is.  In  view  of  this  behavior,  a  new  structure  including  an  additive  material 
with  high  thermal  conductivity  is  designed  for  recovering  cost-free  energy  resources  with  ATEC. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

As  a  result  of  energy  crisis  and  environmental  pollution  fol¬ 
lowed  by  immoderately  burning  fossil  fuels  in  recent  several  dec¬ 
ades,  many  efforts  have  been  poured  into  the  researches  of 
renewable  clean  energy  resources.  Thermoelectric  power  genera¬ 
tion,  because  of  its  simplicity,  less  maintenance  cost,  high  reliabil¬ 
ity  and  environmental  friendliness,  has  become  a  promising  clean 
energy  conversion  technology.  However,  due  to  the  relatively  low 
thermal-to-electric  efficiency  compared  to  conventional  power 
generation  technologies  [1  ],  the  commercialization  of  this  technol¬ 
ogy  is  limited,  mainly  in  the  fields  of  low  cost  or  cost-free  energy 
resources,  such  as  solar  energy  [2-4]  and  waste  heat  from  automo¬ 
tive  exhaust  gases  [5,6], 

The  performance  of  thermoelectric  generator  (TEG)  composed 
of  a  series  of  thermoelectric  couples  is  directly  related  to  the 
dimensionless  figure  of  merit  (ZT),  geometric  arrangement  of 
thermocouple,  external  load  and  temperatures  of  cold  and  hot 
junctions.  To  increase  the  competitiveness  of  thermoelectric 
power  generation  technology,  extensive  investigations  have  been 
conducted. 
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Considerable  improvement  about  ZT  was  obtained  in  recent 
decades.  Zhao  et  al.  [7]  reported  a  ZT  of  2.9  at  923  K  in  SnSe  single 
crystals,  which  was  the  highest  up  to  present.  By  introducing  shape 
parameter,  Sahin  and  Yilbas  [8]  and  Ali  et  al.  [9]  investigated  the 
effect  of  leg  geometry  on  the  performance  of  thermoelectric  cou¬ 
ple.  Ref.  [8]  showed  that  the  efficiency  of  thermoelectric  couple 
increased  notably  when  the  leg  geometry  became  a  trapezium 
shape  along  the  leg  height.  Optimal  external  load  for  peak  power 
of  a  TEG  in  solar  application  was  carried  out  by  Lesage  et  al.  [10] 
using  two  experimental  apparatuses.  They  indicated  that  the  load 
matching  condition  cannot  yield  peak  power  for  any  of  the  thermal 
conditions  tested.  Also,  a  simplified  normalized  power  relation  was 
put  forward. 

The  temperatures  of  cold  and  hot  junctions  were  determined  by 
heat  transfer  capabilities  between  TEG  and  heat  source  or  heat 
sink.  Some  appropriate  approaches  for  increasing  heat  transfer 
capabilities  were  proposed  and  investigated.  In  detail,  Crane  and 
Jackson  [11]  created  a  numerical  model  and  studied  an  integrated 
thermoelectric  cross  flow  heat  exchanger.  Based  on  theoretical  and 
experimental  studies,  some  suggestions  for  performance  improve¬ 
ment  of  TEG  were  proposed  by  Gou  et  al.  [12],  such  as  expanding 
heat  sink  surface  area  in  a  proper  range  and  enhancing  cold-side 
heat  transfer  capacity  in  a  proper  range.  The  performance  of  TEG 
at  various  operating  conditions  was  studied  experimentally  by 
Chen  et  al.  [13],  Results  showed  that  the  effects  of  flow  pattern 
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Nomenclature 

A 

cross-section  area  of  leg,  m2 

/ 

electrical  current,  A 

Greek  symbols 

k 

thermal  conductance  of  leg,  W/(m  K) 

a 

Seebeck  coefficient,  V/K 

K 

overall  thermal  conductance  of  ATEC,  W/(m  K) 

X 

thermal  conductivity  of  leg,  W/(m  K) 

P 

output  power  of  ATEC,  W 

a 

electrical  conductivity  of  leg,  Q  1  m  1 

Q 

heat  flow  rate,  W 

5 9 

angle  in  the  direction  of  circumference,  rad 

Qin 

Fourier  heat  input,  W 

thickness  of  leg,  m 

Qout 

Fourier  heat  output,  W 

6 

dimensionless  temperature  ratio 

Qj 

Joule  heat,  W 
radial  direction 

m 

efficiency  of  ATEC 

R 

Rl 

total  electrical  resistance  of  ATEC,  C2 

external  load,  C2 

annular  shaped  parameter 

Subscripts 

n  n-type  leg 

d  D-tvDe  les 

T 

temperature,  K 

0 

reference  position 

Tm 

average  temperature,  K 

hot  junction 
cold  junction 

ZT m 

the  dimensionless  figure  of  merit 

2 

of  heat  sink  and  water  flow  rate  were  not  significant,  but  the  heat 
source  temperature  played  an  important  role.  Love  et  al.  [14] 
investigated  the  influences  of  heat  exchanger  material  and  fouling 
on  TEG  exhaust  heat  recovery.  They  pointed  out  that  the  perfor¬ 
mance  of  thermoelectric  device  with  stainless  steel  heat  exchanger 
was  lower  than  that  with  aluminum  heat  exchanger  due  to  the  dif¬ 
ference  of  thermal  conductivity,  and  heat  exchanger  fouling  can 
degrade  the  performance  of  studied  devices.  Using  exhaust  gas  of 
vehicles  as  heat  source,  a  mathematical  model  for  the  performance 
of  TEG  devices  was  developed  by  Wang  et  al.  [15],  The  findings 
revealed  that  the  benefit  caused  by  changing  the  convection  heat 
transfer  coefficient  of  high-temperature  side  was  more  pro¬ 
nounced  than  that  of  low-temperature  side.  The  power  output 
characteristics  of  TEG  with  plate-fin  heat  absorber  were  investi¬ 
gated  by  Jang  et  al.  [16]  using  three-dimensional  numerical 
method.  Results  indicated  that  the  optimal  fin  height  and  flue 
gas  velocity  existed  for  maximizing  net  power  density.  In  the 
experiment  conducted  by  Lesage  et  al.  [17],  the  effects  of  two  flow 
channel  inserts  (panel  and  spiral  inserts)  strengthening  the  heat 
transfer  coefficient  on  power  generation  of  liquid-to-liquid  TEG 
were  considered.  It  was  shown  that  the  panel  inserts  enhanced 
power  up  to  110%,  while  the  spiral  inserts  produced  negligible 
power  enhancement.  Using  the  same  experimental  apparatus  pre¬ 
sented  in  Ref.  [17],  the  net  power  output  of  TEG  under  various 
panel  densities  of  flow  impeding  geometric  inserts  was  explored 
by  Amaral  et  al.  [18],  Optimal  insert  panel  density  maximizing 
the  net  power  was  identified  with  respect  to  flow  rate.  A  micro- 
TEG  was  proposed  and  its  performance  was  researched  by  Wojtas 
et  al.  [19],  This  device  had  the  advantages  of  high  heat  transfer 
coefficient  and  small  pumping  power  in  a  compact  volume.  Results 
suggested  that  the  optimization  for  thermoelectric  material  should 
focus  more  on  the  power  factor  than  on  ZT  itself  when  systems  had 
good  thermal  coupling.  Besides,  Wang  et  al.  [20]  designed  the 
geometry  of  heat  sink  such  as  fin  spacing  and  length  with 
two-stage  optimization.  Jang  and  Tsai  [21  ]  optimized  TEG  modules 
pacing  and  its  spreader  thickness  using  a  simplified  conjugate- 
gradient  method.  Tzeng  et  al.  [22]  discussed  the  configurations 
and  dimensions  of  heat  absorber  and  heat  sink  by  applying 
one-dimensional  steady  model.  Also,  heat  pipe  [23,24]  and 
thermosyphon  [25]  heat  sink  were  employed  in  solar  hybrid 
systems  with  TEGs  to  enhance  heat  transfer  capability. 

From  Refs.  [11-25],  two  features  have  been  found.  On  one  hand, 
all  of  them  were  limited  to  flat-plate  thermoelectric  couple  (FTEC), 
where  the  cross-section  shape  of  legs  was  flat-plate.  On  the  other 


hand,  almost  all  of  these  studies  were  concentrated  on  methods 
to  increase  heat  transfer  coefficient,  which  included  altering  flow 
patterns  [13]  and  heat  exchanger  materials  [14],  adding  flow  chan¬ 
nel  inserts  [17,18],  applying  micro  fluidic  heat  transfer  system  [19], 
heat  pipe  [23,24]  and  thermosyphon  [25],  However,  in  some  prac¬ 
tical  applications,  such  as  converting  heat  from  automotive 
exhaust  gases  [5,26]  and  from  coal-fired  boiler  [27],  utilizing  solar 
energy  with  heat  pipe  [23,24]  or  thermosyphon  [25],  the  heat 
source  or  heat  sink  is  cylindrical  in  shape.  In  these  situations,  if 
FTEC  was  still  adopted,  poor  heat  transfer  capability  caused  by 
the  relative  geometries  of  heat  source  or  heat  sink  and  thermoelec¬ 
tric  couple  would  be  presented  [10]  in  spite  of  taking  above  meth¬ 
ods  to  increase  heat  transfer  coefficient.  Consequently,  the 
performance  of  thermoelectric  couple  would  be  degraded  and 
the  commercialization  of  thermoelectric  power  generation  tech¬ 
nology  would  be  further  postponed.  In  this  regard,  annular  ther¬ 
moelectric  couple  (ATEC)  with  annular  shaped  leg  is  proposed, 
ascribing  to  the  fact  that  ATEC  can  contact  with  heat  source  or  heat 
sink  closely  showing  much  smaller  heat  transfer  resistance  than 
that  of  FTEC.  Currently,  the  annular  shaped  thermoelectric  module 
has  already  been  fabricated  [28],  Yet,  to  the  best  knowledge  of 
authors,  the  study  related  to  ATEC  has  received  very  little  atten¬ 
tion,  and  the  impact  of  annular  shaped  parameter  characterizing 
the  feature  of  ATEC  on  the  performance  of  thermoelectric  couple 
is  blank.  What’s  more,  considering  the  fact  that  non-constant 
cross-section  along  thermocouple  leg  and  annular  shape  are  pre¬ 
sented  in  ATEC,  some  new  behaviors  may  occur,  such  as  the  alloca¬ 
tion  of  generated  Joule  heat  on  both  junctions  of  ATEC.  Thus  the 
previous  formulas  preparing  for  FTEC  may  be  inapplicable  to  ATEC. 

In  the  present  study,  taking  ATEC  as  object,  the  fundamental  for¬ 
mulas  of  ATEC  including  temperature  distribution  along  the  legs, 
heat  absorbed  at  hot  junction  and  dissipated  at  cold  junction,  out¬ 
put  power  and  efficiency,  are  derived  based  on  one-dimensional 
steady  model.  Then  the  effect  of  annular  shaped  parameter  on  the 
performance  of  ATEC  under  different  operating  temperature  ratios, 
external  loads  and  ZT  is  theoretically  analyzed.  In  the  end,  some 
suggestions  or  notes  for  practical  applications  of  ATEC  are  provided. 

2.  Theoretical  model 

A  typical  ATEC  is  presented  in  Fig.  1,  which  consists  of  a  p-  and 
n-type  thermoelectric  legs  sandwiched  between  heat  source  and 
heat  sink.  The  first  part  of  this  study  is  to  derive  the  fundamental 
formulas  of  ATEC.  Fig.  2  shows  the  physical  model  of  one  leg.  In 
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M(r) 


dA(r)  dT  I2 
dr2  T  dr  dr  +  oA(r)  ~ 


(2) 


where  I  is  the  current  following  through  leg.  2,  a  and  T  are  the  ther¬ 
mal  conductivity,  electrical  conductivity  and  temperature  of  leg, 
respectively.  A(r)  is  the  cross-section  area  of  leg,  which  can  be 
expressed  as 


A(r)  =  Acp  ■  r  ■  5 


(3) 


Taking  the  area  A(r0)  in  the  position  of  r0  =  (r,  +  r2)/ 2  as  refer¬ 
ence,  A{r)  can  be  expressed  as 


A(r)=A(r0) 


to 


(4) 


where  A(p  and  <5  are  respectively  the  angle  in  the  direction  of  cir¬ 
cumference  and  thickness  of  leg.  ri  and  r2  represent  the  positions 
of  hot  and  cold  junctions  of  ATEC.  Substituting  Eq.  (4)  into  Eq.  (2), 
Eq.  (2)  can  be  reformed  as 


dr2  r  dr  KrJ  la A2(r0) 

According  to  the  assumptions,  the  boundary  conditions  are 

T(r=rr)  =  Ti,  T(r  =  r2)  =  T2  (6) 

To  reduce  the  cost  of  manufacture,  the  geometric  dimensions  of 
p-  and  n-type  legs  are  the  same  as  Ref.  [29],  i.e.,  Ap(r)  =An(r)  =A(r), 
A<pp  =  A <pn  =  Acp,  5P  =  Sn  =  <5,  rlp  =  rln  =  r,,  r2p  =  r2n  =  r2.  Hence,  the 
temperature  distributions  of  p-  and  n-type  legs  are  similar  and 
can  be  expressed  respectively  as 


Tp(r) 


I2rl 


.  (Inr) 


22perpA2(r0) 

|72r2(lnri+lnr2)  T,  -T2  1  ] 

[  22p(j pA2(r0)  +Bri-tor2J 

T2lnri  -  Ti  lnr2  f2rglnr,lnr2 
lnr,  -  lnr2  22po-pA2(r0) 


J2r2 


-(lnr)2 


the  theoretical  analysis,  the  following  approximations  are  made 
[2,8,9]: 

•  The  properties  of  leg  material  are  independent  to  temperature 
and  position,  so  the  Thomson  effect  is  neglected. 

•  The  temperatures  of  hot  and  cold  junctions  are  maintained  at 
constant  values  of  T,  and  T2  respectively. 

•  Electric  and  thermal  contact  resistances  are  negligible. 

•  One-dimensional  heat  transfer  along  the  radial  direction  r  is 
considered,  and  the  lateral  surfaces  are  considered  to  be 
adiabatic. 

•  Along  the  leg  (p-  or  n-type  leg)  height  direction,  the  leg  has  the 
same  thickness. 

These  assumptions  allow  us  to  focus  on  the  fundamental  infor¬ 
mation  of  ATEC.  Applying  energy  balance  to  the  infinitesimal  ele¬ 
ment  dr  shown  in  Fig.  2,  it  has 

Qin  +  Qj  -  Qout  =  0  (1) 

where  Qin,  Qout  and  Qj  are  the  Fourier  heat  input,  Fourier  heat  output 
and  generated  Joule  heat,  respectively.  Eq.  ( 1 )  can  be  rewritten  as 


22„crnA  (r0) 

rga(lnr,+lnr2)  t,-r,  1 

[  2/nfjnA2(r0)  lnr,  -  lnr2J 
T2  lnr,  -  T,  lnr2  J*rglnrt  lhr2 
lnr,  -  lnr2  22n(7 nA2(r0) 

The  heat  transferred  through  the  leg  of  ATEC  along  r  is  pre¬ 
sented  as 


(,0) 

Thus  the  thermal  conductance  k  of  leg  can  be  defined  as 

k=  ;'A<P.5  (11) 

In  r2  -  In  r,  v  ’ 

The  overall  thermal  conductance  K  of  ATEC  shown  in  Fig.  1  is 

K=kp  +  kn  =  +  (12) 

F  In  r2  -  In  X\  In  r2  -  In  n  v  x 
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On  the  other  hand,  the  total  electrical  resistance  R  of  ATEC  can 
be  presented  as 

R=  r  (  dr  dr  ^  =  lnr2  - Ino  lnr2  - lnr, 

Jr,  \ffpA(x)  MWj  ov  \(p-5  an  A(p  S  [  ’ 
The  heat  absorbed  at  the  hot  junction  Q,  and  dissipated  at  the 


cold  junction  Ch  of  ATEC  are 

Qi=(ap-a  n)JT1-[^A(r)^]r 

(14) 

Q2  =  (ap-an)/T2-[2pA(r)^]r 

~hr)f L 

(15) 

where  ap  and  an  are  the  Seebeck  coefficients  of  p-  and  n-type  legs 
respectively. 

Now  combining  Eqs.  (4),  (7),  (8),  (12),  and  (13)  with  Eqs.  (14) 
and  (15),  Q,  and  Q2  become 


Q,  =  a  IT,  +  K(T !  -  T2)  -  0.5 12R 
Q_2  =  a IT2  +  I((T  1  -  T2 )  +  0.5 12R 


(16) 

(17) 


where  a  =  ap  —  an  is  the  Seebeck  coefficient  of  ATEC.  When  the 
external  resistance  RL  is  imposed  on  ATEC,  the  current  I  through  it 
can  be  expressed  as  follows 


_  X„  ■  Acp  •  5  •  (r,  +r2)/ 2 
r2  -  r  . 


(23) 

(24) 

(25) 


"  <T„  ■  A(p  •  5  •  (rj  +  r2)/ 2 
The  annular  shaped  parameter  sr  is  defined  as 
Sr  =  r2/r, 

Further  assuming  the  similar  doped  alloys  are  used  to  make  p- 
and  n-type  ATEC  legs,  so  Xp  =  Xn  =  X,  erp  =  <xn  =  a,  ap  =  -an  =  0.5a  [29], 
Combining  this  hypothesis  with  Eqs.  (12),  (13),  and  (23),  (24),  (25), 
the  dimensionless  forms  of  R  and  K  can  be  obtained  as 


(26) 


R  _  (sr  +  1)  lnsr 
Ro~  sr- 1 

M-V  (27) 

K0  (Sr  +  f}4tasr  {  > 

Substituting  Eqs.  (21  )-(24),  (26),  and  (27)  into  Eqs.  (19)  and  (20), 
the  dimensionless  forms  of  P  and  r\,  which  are  the  implicit  functions 
of  sT,  can  be  written  as 

P  2(K/Ko)(ZTm)(RL/R0)(R/R0)(Q  -  l)2 

w2  (0  +  i)[(rl/ro)  +  (R/ro)]2 


2(ZTm)(0-l)(RL/Ro)(R/Ro) 


(' e  +  1)[(Rl/R0)  +  (R/R0))2  +  20(ZTm)(i?/Ro)[(RL/Ro)  +  (R/Ro)]  -  ( ZTm)(0  -  1  )(R/R0f 


«(Ti-T2) 

-  Rib'  (18 

The  output  power  P  and  efficiency  rj  of  an  ATEC  can  be  given  a: 


Equating  the  first  partial  derivative  of  Eq.  (28)  with  respect  to  RL 
to  zero,  RL  corresponding  to  the  maximum  output  power  can  be 
obtained  as  follows 


=  _ «2(Ti  —  T2)Rl _ 

0.1  a2T,(R  +  RL)  +  K(R  +  Rl)2  -  0.5a 2(T,  -  T2)R 


(20) 


As  can  be  seen,  the  forms  of  Eqs.  (16),  (17),  (19),  and  (20)  are  the 
same  as  those  of  FTEC  with  constant  cross-section  area  [30],  And 
from  the  above  derivation  process,  it  can  be  concluded  that  for 
any  types  of  thermoelectric  couple  (FTEC,  ATEC  or  others)  and 
cross-section  geometries,  the  corresponding  forms  of  Qj,  Q2,  P 
and  r\  must  be  the  same  except  for  the  expressions  of  K  and  R, 
provided  that  the  mentioned  assumptions  are  satisfied.  It  provides 
the  theoretical  basis  for  Refs.  [8,9],  where  the  source  of  formula 
about  the  efficiency  of  FTEC  with  non-constant  cross-section  is 
not  given. 

The  dimensionless  figure  of  merit  ZTm  based  on  the  average 
temperature  Tm  =  (T,  +  T2)/2  is  provided  as 


ZTm 


a2  (T1+T2)_a2(0+1)T2 
~  KR~  2  2  KR 


Eq.  (30)  agrees  with  Ref.  [2],  Similarly,  RL  for  maximum  effi¬ 
ciency,  which  is  the  same  as  that  in  Ref.  [2],  is  given  by 


Then  the  maximum  dimensionless  output  power  and  efficiency 
of  ATEC,  respectively,  are 


(&L- 


(ZTm)(K/K0)(e-\)2 

2(0+1) 


2(Zrm)(fl-l)vT+Z7^ 


max  (0+l)(l+v/T+Z7m)2  +  20(ZTm)(l  +  v/T+ZTi;)-(Zrm)(0-l) 

(33) 

According  to  Ref.  [31],  Eq.  (33)  can  be  further  reformed  as 


(21) 


^7 max  =  - 


0-1  vT+ZTm  - 1 


(34) 


where  0  is  the  dimensionless  temperature  ratio  and  defined  as 

0  =  r,/r2  (22) 

The  referenced  thermal  conductance  K0  and  electrical  resistance 
R0  based  on  the  position  r0  of  n-type  leg  [8],  respectively,  can  be 
expressed  as 


Eqs.  (32)-(34)  clearly  show  that  in  ATEC,  the  maximum  dimen¬ 
sionless  output  power  is  affected  by  sr,  whereas  the  maximum  effi¬ 
ciency  remains  unchanged  with  varying  sr.  Furthermore,  the 
conclusion  that  the  geometries  of  thermoelectric  couple  legs  have 
no  effect  on  the  maximum  efficiency  can  be  inferred  according  to 
Eqs.  (33)  and  (34). 
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Fig.  3.  Variation  of  P/(K0r2)  with  sr  at  different  9  for  RJRq  =  4.0,  ZTm  =  1.5. 


It  should  be  noted  that  sr  <  1  represents  a  kind  of  ATEC  that  the 
cold  junction  locates  in  the  concave  side,  while  the  hot  junction  is 
placed  in  the  concave  side  for  sr  >  1  (like  the  situation  shown  in 
Fig.  1).  ATEC  is  transformed  into  FTEC  on  the  condition  that 
sr  =  1,  i.e.,  both  r,  and  r2  tend  to  be  infinite.  In  this  case,  the  derived 
formulas,  like  Eqs.  (16)  and  (17),  are  the  same  as  Ref.[30]  totally, 
demonstrating  that  the  above  theoretical  model  is  valid. 

Due  to  the  fact  that  the  leg  length  of  ATEC  compared  to  ri  or  r2  is 
relatively  small,  the  considered  range  of  sr  is  between  0.1  and  4.0. 
Meantime,  the  following  parameters  are  employed:  0  =  2.0, 
RJR0  =  4.0  and  ZTm  =  1.5. 

3.  Results  and  discussion 

Under  the  conditions  of  different  0,  R\J Ro  and  ZTm,  theoretical 
analysis  of  ATEC  is  conducted  to  examine  the  effect  of  sr  on 
PKK0T2)  and  rj  of  ATEC.  The  leg  length  (r2-r,)  and  A cp  are  assumed 
to  be  constant. 

Fig.  3  shows  variation  of  P!{K0T2)  with  sr  for  different  0.  There 
exists  a  maximum  for  P/(ff0T2)  as  sr  varies,  and  the  corresponding 
sr  is  equal  to  1.  This  suggests  that  the  performance  of  FTEC  is  more 
excellent  than  that  of  ATEC  from  the  perspective  of  PI{K0T2).  It  is 
expected  and  the  reasons  are  presented  as  follows.  For  sr<l, 
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although  much  heat  can  be  absorbed  at  the  hot  junction,  the  rela¬ 
tively  small  heat  transfer  area  at  the  cold  junction  limits  PI{KqT2)  of 
ATEC.  And  in  the  case  of  sr  >  1 ,  since  the  absorbed  heat  at  the  hot 
junction  is  diminished  compared  to  FTEC,  PI(K0T2 )  is  also  less  than 
that  of  sr=  1.  This  behavior  can  also  be  explained  as:  in  the 
condition  of  constant  Seebeck  voltage  (resulting  from  the  constant 
temperatures  of  junctions  Ti  and  T2),  as  the  farther  sr  is  away 
from  1,  the  larger  the  electric  resistance  (Eq.  (26))  is,  hence,  the 
lower  PI(K0T2)  is.  However,  in  some  applications  of  heat  source 
or  heat  sink  with  round  shape  [5,23-27],  heat  transfer  capability 
between  thermoelectric  couple  and  heat  source  or  heat  sink  is 
reduced  considerably  by  geometries  mismatching,  provided  that 
thermoelectric  couple  with  suitable  leg  shape  is  not  employed. 
Consequently,  the  performance  of  thermoelectric  couple  may  be 
deteriorated  seriously.  This  fact  indicates  that  FTEC  may  not  be 
the  best  choice  in  these  applications  and  the  selection  of  thermo¬ 
electric  couple  with  proper  leg  shape  is  of  great  importance. 
Fig.  3  also  shows  that  increasing  0  does  not  affect  the  variation  of 
PI(KqT2)  with  sr.  For  example,  as  sr  varies  from  1.0  to  2.0,  PI{K0T2 ) 
is  dropped  by  2.6%  for  all  considered  0. 

Variation  of  rj  with  sr  for  different  0  is  shown  in  Fig.  4.  sr  approx¬ 
imately  has  no  effect  on  r\  in  the  studied  case  of  RL/R0  =  4.0  and 
ZTm  =  1.5.  This  is  attributed  to  the  values  of  considered  RL/Rn  and 
those  of  Rl/Ro  for  maximizing  rj.  In  detail,  since  the  center  section 
r0  of  n-type  leg  is  considered  as  the  dimensionless  benchmark,  R/ 
R0  has  approximately  the  value  of  2.0.  According  to  Eq.  (31),  RL/ 
R0  for  maximizing  t]  equals  approximately  to  3.2,  which  has  no 
considerable  difference  with  the  considered  Rl/Ro  (Rl/Ro  =  4.0). 
Therefore,  in  this  situation,  as  shown  in  Eqs.  (33)  and  (34),  the 
influence  of  sr  on  rj  is  nearly  negligible.  This  behavior  between  sr 
and  rj  disappears  when  RL/Ro  deviates  from  the  condition  for  max¬ 
imizing  >y,  which  can  be  seen  from  Fig.  6. 

Fig.  5  presents  the  variation  of  PI{K0T2)  with  sr  under  different 
Rl/Ro-  A  maximum  of  PI{K0T2 )  is  detected  when  equaling  sr  to  1. 
Namely,  FTEC  outperforms  than  other  typed  thermoelectric 
couples  with  sr  #  1.  This  behavior  is  smeared  out  with  increasing 
Rl/Rq.  As  the  result  presented  in  Eq.  (30),  Fig.  5  shows  that  the  max¬ 
imum  of  PI(K0T2 )  occurs  at  Rl/Ro  =  2.0.  The  effect  of  sr  on  >;  under 
different  RJR0  is  given  in  Fig.  6,  where  some  interesting  results 
are  found.  Specifically,  with  RL/Ro  =  2.0,  as  sr  increases,  r]  increases 
before  attaining  a  maximum  after  which  it  decreases,  while 
this  behavior  reverses  and  a  minimum  of  rj  emerges  when 
Rl/Ro  5s  4.0.  The  corresponding  sr  for  these  extremums  are  sr=  1. 
Besides,  as  long  as  RL  satisfies  the  condition  shown  in  Eq.  (31),  t] 


Fig.  4.  Variation  of  rj  wit 
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Fig.  7.  Variation  of  P/(KoT2)  with  sr  at  different  ZTm  for  RJRn  =  4.0,  B  =  2.0. 


sink.  However,  in  applications  of  heat  source  or  heat  sink 
with  round  shape  [5,23-27],  thermoelectric  couple  with 
suitable  leg  shape,  i.e.,  ATEC,  is  obligatory  to  eliminate  the 
heat  transfer  resistance  caused  by  the  geometries  mismatch¬ 
ing.  Therefore,  to  improve  the  performance  of  thermoelectric 
couple,  the  effective  contacting  area  (for  instance,  using 
thermoelectric  couple  in  an  appropriate  shape)  should  be 
increased  simultaneously  with  the  increment  of  heat  trans¬ 
fer  coefficient. 

(3)  Figs.  3, 5  and  7  show  that  the  dimensionless  power  output  of 
ATEC  with  sr  =  1  (FTEC)  is  higher  than  those  with  sr  #  1.  Yet, 
previous  analysis  demonstrates  that  ATEC  must  be  used  for 
round  shaped  heat  source  or  heat  sink.  Combining  these 
two  points  as  far  as  possible,  a  new  structure  presented  in 
Fig.  9  is  proposed  for  recovering  cost-free  energy  resources, 
such  as  exhaust  waste  heat  and  solar  energy.  In  this  struc¬ 
ture,  an  additive  material  with  high  thermal  conductivity 
is  employed  between  ATEC  and  round  shaped  heat  source 
or  heat  sink.  The  advantages  of  this  structure  are  listed  as 
follows.  Firstly,  compared  to  previous  structure  without 


will  reach  the  maximum  regardless  of  sr.  These  results  imply  that 
for  practical  applications  of  ATEC,  much  attention  should  be  paid 
to  the  connected  external  load. 

Figs.  7  and  8  show  the  variations  of  PI(K0T2)  and  r\  with  sr  for  dif¬ 
ferent  ZTm.  Since  the  behaviors  between  PI(KqT2)  and  sr  as  well  as 
between  r]  and  sr  are  similar  to  those  under  different  0  (shown  in 
Figs.  3  and  4),  the  descriptions  are  not  given  in  detail  for  the  pur¬ 
poses  of  brevity. 

To  enhance  the  practicability  of  ATEC  and  to  facilitate  the  prac¬ 
tical  applications,  some  suggestions  or  notes  are  provided. 

(1)  Above  analysis  indicates  that  compared  to  8  and  ZTm,  the 
influence  of  sr  on  the  performance  of  ATEC  is  affected  by 
RiJRo  seriously.  Hence,  enough  attention  should  be  paid  to 
the  external  load  matching  when  applying  ATEC. 

(2)  Heat  transfer  capability  between  thermoelectric  couple  and 
heat  source  or  heat  sink,  as  a  critical  parameter  for  the  per¬ 
formance  of  thermoelectric  couple,  is  determined  by  both 
heat  transfer  coefficient  and  corresponding  effective  con¬ 
tacting  area.  Extensive  work  has  been  done  mainly  to 
enhance  the  heat  transfer  coefficient,  while  for  effective  con¬ 
tacting  area,  only  little  work  has  been  done.  Such  a  situation 
is  probably  understandable  for  flat-plate  heat  source  or  heat 


Round  shaped  heat  reservoir 


Fig.  9.  Schematic  of  the  proposed  structure  (only  one  ATEC  is  presented). 
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the  additive  material,  sr  in  this  structure  is  more  close  to  1, 
and  consequently  ATEC  in  this  structure  has  better  perfor¬ 
mance.  Also,  the  thermal  stress  is  reduced,  which  can  extend 
the  lifetime  of  ATEC.  Secondly,  the  effective  contacting  area, 
then  heat  transfer  capability  is  increased  compared  to  using 
FTEC.  Consequently  the  performance  of  ATEC  is  improved. 
However,  since  an  extra  contacting  surface  between  ATEC 
and  round  shaped  heat  source  or  heat  sink  is  added  by 
comparison  to  the  previous  structure  without  the  additive 
material,  care  should  be  taken  into  account  when  using  this 
structure. 


4.  Conclusions 

To  enhance  heat  transfer  capability  from  the  aspect  of  increas¬ 
ing  effective  contacting  area,  the  annular  thermoelectric  couple 
(ATEC)  is  proposed  for  the  practical  applications  with  round 
shaped  heat  source  or  heat  sink.  An  annular  shaped  parameter  sr 
is  utilized  to  characterize  ATEC.  The  fundamental  formulas  about 
ATEC  are  derived  based  on  one-dimensional  steady  model. 
Theoretical  analysis  of  the  effect  of  sr  on  the  performance  of  ATEC 
under  different  temperature  ratios,  external  loads  and  dimension¬ 
less  figure  of  merits  is  carried  out. 

The  results  show  that  the  fundamental  formulas  of  ATEC  are 
similar  to  those  of  flat-plate  thermoelectric  couple  (FTEC)  except 
for  the  expressions  of  overall  thermal  conductance  and  total  elec¬ 
trical  resistance.  Specially,  as  sr  equals  to  1 ,  ATEC  turns  into  FTEC. 
The  extremums  of  dimensionless  output  power  and  efficiency  are 
obtained  when  sr  =  1.  The  influence  of  sr  on  the  dimensionless  out¬ 
put  power  and  efficiency  of  ATEC  is  affected  by  the  external  load 
seriously  as  compared  to  the  temperature  ratio  and  dimensionless 
figure  of  merit.  As  external  load  increases,  the  extremum  type  of 
efficiency  shifts  from  maximum  to  minimum.  While  for  dimen¬ 
sionless  output  power,  the  type  is  always  the  maximum,  but  the 
effect  of  sr  becomes  negligible.  These  imply  that  the  selection  of 
the  external  load  for  ATEC  is  of  great  importance.  For  generating 
power  with  cost-free  energy  resources,  considering  the  fact  that 
ATEC  with  sr  =  1  has  higher  dimensionless  output  power  than  those 
with  sr  #  1,  a  new  structure  for  ATEC  is  designed.  In  this  structure, 
an  additive  material  with  high  thermal  conductivity  between  ATEC 
and  round  shaped  heat  source  or  heat  sink  is  employed. 

A  more  accurate  theoretical  model  will  be  developed  in  the  near 
future,  in  which  the  temperature-dependant  properties  of  thermo¬ 
electric  material,  Thomson  effect  and  non-adiabatic  condition  on 
lateral  surface  of  ATEC  legs  will  be  considered. 
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